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Q I Abstract 

f^ ■ We analyze dynamical consequences of a conjecture that there ex- 

ists a fundamental (indivisible) quant of time. In particular we study 
0^] the problem of discrete energy levels of hydrogen atom. We are able 

to reconstruct potential which in discrete time formalism leads to en- 
ergy levels of unperturbed hydrogen atom. We also consider linear 



^ , energy levels of quantum harmonic oscillator and show how they are 

^^' produced in the discrete time formalism. More generally, we show that 

in discrete time formalism finite motion in central potential leads to 
K> , discrete energy spectrum, the property which is common for quantum 

\—i I mechanical theory. Thus deterministic (but discrete time!) dynamics 

is compatible with discrete energy levels. 

1 Introduction 

Discovery of discrete energy levels for atoms demonstrated that the classical 
Newtonian model could not be used to describe this phenomenon. Now 
days this phenomenon is described in the framework of quantum mechanics. 
The advantages of the computational methods of quantum mechanics are 
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well known. They are widely used for computation of energy levels not 
only of atoms, but also in essentially more complicated situations. However, 
despite this computational power, quantum mechanics induced many still 
unsolvable phenomenological problems, see PO-[SZI for extended discussions. 
One of distinguishing features of quantum mechanics (at least in Copenhagen 
interpretation) is the impossibility to provide realist deterministic description 
of quantum reality. In particular, here particles do not have well defined 
trajectories. 

The impossibility to use the deterministic evolutionary model looks rather 
counterintuitive. Moreover, it contradicts to many quantum experiments 
where "trajectories" of particles are well observed, e.g. in Wilson's camera. 
Of course, in modern quantum phenomenology this problem is solved via 
the principle of complementarity. Nevertheless, this strong deviation from 
our intuitive picture of physical reality induces a rather general opinion that 
quantum mechanics has a lot of mysteries |27] . 

This unsatisfactory status of quantum mechanics induces new and new 
attempts to create new models that would be closer to our physical intuition. 
We just mention one model, Bohmian mechanics pm i21j. 

Last years there were performed intensive investigations to reconsider 
probabilistic foundations of quantum mechanics, see e.g. [321 EII- These 
investigations are of the great importance. It is well known that "quantum 
probability" differs strongly from "classical probability". Typically, see e.g. 
^27\ . it is pointed out that quantum randomness is irreducible and funda- 
mental (in the opposite to classical randomness that could be reduced to e.g. 
randomness of initial conditions). This irreducible quantum randomness is 
deeply connected with the impossibility of realist deterministic models that 
would reproduce quantum probabilities. It seems to be impossible to imagine 
that quantum particles have deterministic trajectories, since the existence of 
trajectories should imply classical probabilistic rules. But these rules are 
violated (for example, in the two-slit experiment J2Z|; see also investigations 
on EPR-Bohm-Bell consideration and chameleon effect PTIITHI IT^). 

In papers j^Hl ISHl HOI HI] there was developed the contextual approach 
to quantum probabilities that might be used to explain the origin of quan- 
tum randomness in the classical (but contextual!) probabilistic framework. 
Therefore the probabilistic constraint to create realist deterministic models 
need not be taken into account^. 



^Of course, there are also Bell's inequality arguments. But there arguments are strongly 



What kind of realist deterministic models could be created for quantum 
phenomena? 

As we have already mentioned, we could not directly use Newtonian me- 
chanics even for the simplest quantum system - the hydrogen atom. Thus 
Newton equations should be changed. The most straightforward idea is to 
change classical forces (in the case of atom Coulomb's law) and find new 
forces, quantum forces. Such forces should "drive" particles along trajecto- 
ries that reproduce quantum data. There is a rather common view point 
that one of successful realizations of such a program is given by Bohmian 
mechanics. However, this problem is not simple, since Bohmian mechanics 
is, in fact, not mechanics, but a field model. There exists an additional field 
equation for quantum potential. We do not know any model that would give 
the realization of the discussed program of "modernization" of Newtonian 
mechanics. 

In the present paper we consider discrete time Newtonian model. This is 
a kind of classical physical model (in particular, realist deterministic). The 
only difference is discreteness of time. Thus we use the classical force - in- 
teraction picture, but the discrete time version of Newton's second law. One 
of the main advantages of this model is its simplicity. We do not change 
phenomenology of classical physics (position, velocity, force-interaction, tra- 
jectories). We only change the mathematical representation of time. 

On the other hand, discreteness is the main distinguishing feature of quan- 
tum physics. In fact, M. Planck and A. Einstein obtained Wien's law simply 
by assuming discreteness of energy. However, this discreteness approach was 
not developed further to get a formalism based only on the discreteness pos- 
tulate. Discreteness of quantum observables was reproduced by using an 
advanced mathematical formalism based on the representation of the ob- 
servables in the complex Hilbert space. As we have already mentioned, the 
use of this formalism (despite its computational advantages) induced many 
phenomenological problems. We are trying to modify classical physics by 
starting with one natural postulate: 

TD: "Time is discrete". 

Here, in particular, continuous Newton equations (differential equations) 
are just approximations of real equations of motions, namely, Newton's dif- 



coupled to nonlocality. And this problem is far from problems considered in the present 
paper. 



ference equations^. In the present paper we use the discrete time model for 
the hydrogen atom. The discrete-time postulate implies discrete orbits and 
energy levels. At the same time we have deterministic motion along orbits. 
In the limit we get continuous motion along circular orbit (a kind of Bohr's 
correspondence principle). 

In the present paper we do not try to develop some statistical theory 
(a kind of Born approach). Our investigation has some similarities with 
the original paper of W. Heisenberg (that was not statistical one), ^ |2]. 
We hope that starting with only discrete time postulate we would be able 
to develop simpler formalism to calculate spectra of quantum observables - 
without using noncommutative calculus and without a cardinal change of the 
phenomenology of classical physics. 

It is interesting to point that there is some similarity between our ap- 
proach and G. 't Hooft's approach IHl IHl where a general scheme was 
proposed that maps states of quantum field system to the states of a com- 
pletely deterministic field model. Although in this note we do not consider 
the field-particle duality it would be interesting to study this problem. 

Some statistical consequences of the postulate (TD) were investigated in 
our previous papers j4^ 



2 Discrete Time Dynamics 

In classical mechanics a dynamical function A = A{p, q) (here p and q are 
momenta and coordinates of the system) evolves according to the following 
well known equation |3j 

DtA = {A,H} (1) 

where H = H{p, q) is a Hamiltonian of the system and in the right hand side 
is a Poisson bracket, which could be presented as 

{A,B} = —^-^^ (2) 

dq dp dp dq 

The left hand side of (J]) contains a continuous time derivative 

DtA = —- 

dt 

^We reverse the modern viewpoint on the description of physical reahty. Not (discrete) 
difference equations are used to approximate (continuous) differential equations, but the 



inverse 



As it was mentioned earlier we are interested in construction dynamics with 
discrete time. This is done with the help of discrete derivative which is 
postulated to be 

D\^U = -[A{t + T)-A{t)l 

T 

where r is the discreteness parameter. This parameter is finite and is treated 
in the same way as Plank constant in quantum mechanical formalism. In 
particular if r is small relative to dimensions of the system then classical 
approximation with continuous derivative might work well (although this 
could not be the case all the time in the same sense as there are examples 
when quantum formalism is reasonable even for macroscopic systems, for 
example in superfiuidity) . 

Summarizing, the discrete time dynamical equation is postulated to be 

Di-^A = {A,H}, (3) 

where A{p, q) is a real- valued function of real-valued dynamical variables and 
in the right hand side there is classical Poisson bracket . The equation ^ 
could be solved in the sense that we can write 

A{t + T) = A{t) + T{A,H} (4) 

thus providing the evolution of any dynamical function A = A{p, q). 
Note that in our model the coordinate space is continuous. 



3 Motion in Central Potential 

Here we will study the properties of motion in central potential U = U{r) in 
discrete time formalism. As we will see discreteness of time enriches mechan- 
ics with some new properties which are usually thought as having quantum 
nature. In particular as it will be shown below in discrete time mechanics 
stationary orbits (i.e. finite motion) have discrete energy spectrum. We point 
that the phase space is assumed here to be a continuous real manifold. 

Following the general approach described in the previous section we start 
from the classical Hamiltonian and then write the dynamical equations. In 
polar coordinates (r, ip) the Hamiltonian of the system with mass m in central 
potential f/(r) is given by 

i7(r,p„^,p^) = |^ + -i^ + ^(r), (5) 
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Figure 1: First three trajectories. 



n=3 



where pr and p^ denote momenta corresponding to r and if - radial and angu- 
lar coordinates respectively. Using ^ let us write the dynamical equations. 
We obtain 



r{t + T) =r(t) + r 
Prit + T 
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(6) 
(7) 

(8) 
(9) 



The equation Q corresponds to angular momentum conservation in central 
field - this is a direct analog of the angular momentum conservation law in 
classical (continuous time) dynamics and is a consequence of the fact that 
our Hamiltonian does not depend on Lp, (it is a so called cyclic variable). 

Let us limit ourself to circular stationary periodic orbits. In this case 
we should have r{t + t) = r(t) and thus using (jH)) we see that the radial 
momentum should be zero, 

Prit) = (10) 

From ((Zj) and (|TTH) we obtain the following condition 



PZ 



mr^ 



dU_ 
dr 



(11) 



Let us finally come to the angular coordinate (f. For stable motion the 
following periodicity condition should be satisfied (Fig. Q) 



^{nr) = (^(0) + 27r, 



:i2) 



where n = 1, 2, . . . (note that discreteness of n is a consequence of discreteness 
of time). Using (jHI) and (fT^ we get 

nr—^ = 2n, 

or 

27Tmr'^ , . 

P^ = 13 

nr 

From (fTT|) and (fT^ we get the following equation for the radius of the n-th 

orbit 

Aii'^mr dU ,, ,, 

-^T- = IT 14 

If potential U{r) is known then from equation ()14|1 we can find r = r„. If 
we consider physical potentials, i.e. potentials for which the force, — f/'(r), is 
smooth, negative, and is strictly monotonically decreasing in absolute value 
as r grows, vanishing on infinity, then solution of ()14|) always exists and 
unique. Now upon substituting r„ to the Hamiltonian (j3)) we obtain energy 
levels En. 

For circular periodic orbits the original Hamiltonian (jSJ due to |TU|) and 
fjll|) simplifies to the following form 

H = -r-— + U(r) 
2 dr ^ ' 

Thus the expression for energy E^ of the n-th stable periodic orbit in terms 
of its radius r„ is given by 



F -1 ^ 
" 2''" dr 



+ f/(r„), n = l,2,... (15) 



As we see if potential U{r) allows stationary periodic motion and equation 
(HD) has unique positive solutions r„ then the energy spectrum is discrete. 
This situation is directly analogous to quantum mechanics where for finite 
motion we might expect discrete energy levels. 

4 Energy Levels of Hydrogen Atom 

Our task in this section is to study whether the discrete time dynamics 
in central field described can lead to energy levels of hydrogen atom. We 



treat energy levels as given (measured) quantities and our task is to find the 
corresponding potential. We restrict ourself to the simplest case when the 
atom is unperturbed by external electric or magnetic fields and thus currently 
we do not study splitting of the energy levels (note that in order to observe 
"degenerate" levels in quantum mechanical treatment one needs to somehow 
perturb the system in such a way that the levels split). 

We start from the following energy spectrum for hydrogen atom ^3] 



En = K 



n 



2' 



n 



1,2, 



(16) 



where 7 ~ 13.6 eV is ionization energy of the hydrogen atom. The task is to 
find such U = U{r) that leads to the spectrum (fTHjl . Using (|T^ . (fTH|) . and 
(fTB)) we get 



1 r' 
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.j_| + C(,„) = ^^, (17) 

where the constant ^ is given by (note that ^ depends on discreteness pa- 
rameter r) 

? = ^ (18) 



and n = 1, 2, 



Let us rewrite equation ()17|) in the following form 



f/(r. 



n^ 



'-XK + i 



(19) 



We want to find U{r). The idea is to obtain dependence of r„ = f{n) on 



n and then inverting it substitute n = f 



-1; 



in (fT^ . here / ^ denotes 



function inverse to /. As a result of this procedure we will get rid of explicit 
dependence of the right hand side of p9|) on n, it will depend only on r„. 
Then we interpolate the result for any r ^ 0. The resulting U{r) we check 
by substitution to (ITl|l - (ITH|l . 

Let us proceeding as described above. First, we have to find dependence 
of r„ on n. To do this let us assume that n is continuous and take the 
derivative in n of both sides of (fT^ . This gives 



dU 
dr 

Now, from ()14|) we find 



dvn 
dn 

dU 
dr 



S n 



^rnK 



ri'' 



+ H2 

n^ 



n^ 



(20) 



(21) 




Figure 2: The shape of potential ()25|) which leads to energy levels of the 
hydrogen atom (we put /? = ^ = 7 = !)• 

Substituting (pTjl to (j^ we obtain the following differential equation for r„ 



2enr„r^-er^-27 = 



(22) 



Solving (j^ and taking into consideration only positive solution we obtain 



-27 + e2/3«n 



(23) 



where the constant /? is due to integration. Inverting (|^ we get 

n = e-'P^{2^ + rli) (24) 

Substituting (jUI) to ^^ we obtain 



f/(r„^ 



.4/3? 



47 + 2r2e 
or performing interpolation, i.e. putting r in place of r„ we finally get 

p4/3? 
f/(r) : 



47 + 2r2^ 



(25) 



Now, substituting (j^ to (fT^ - (fT3j) we come to the expected energy levels 
fll6|) of hydrogen atom. 

The form of the potential (J25|l is presented on Fig|2l It is interesting to 
note that unlike Coulomb's potential it is nonsingular at r = 0. 



5 Spectrum of Harmonic Oscillator 

Procedure described in the previous section could be used to obtain potentials 
corresponding to arbitrary energy spectrum. Here we deduce the potential 
which results in linear energy levels of the homogeneous two dimensional 
quantum harmonic oscillator. 

In quantum mechanics homogeneous two dimensional harmonic oscillator 
is a system described by the potential 

U{r) = -mtuV^ (26) 

Solving the Schrodinger equation in potential ()26p results in the following 
energy spectrum 

EA = hiu{A + l), A = 0,1,... (27) 

The simplest way to get this relation is to note that we deal with two uncou- 
pled oscillators (indeed, if x and y are Cartesian coordinates then equation 
(PUj) takes the form U = |ma;^(x^-|-?/^)), each having energy E^ = hu{k + j), 
then the total energy is just the sum of two such terms and we get (jTfj) . 

Let us first rewrite ^I7\i in terms of n = 1,2,..., we have n = A + 1 and 
thus 

En = an, n = 1, 2, . . . , (28) 

where a = fiw. Our task is to find potential U = U{r) which result in 
energy levels (PH|) . Proceeding as in the previous section we assume that n 
is continuous and write the differential equation 

^ = " (29) 

Now, using (jl4|) and (jT^ equation ^I^ reduces to 

/ 2 

ryt ly ' ty'-' 

2e^ - e^ = «, (30) 



where ^ is given by (jTSj). The positive solution of ^U^ is given by 



an^ 



'/?^+^, (31) 
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where /3 is a constant due to integration. This leads us to a cubic equation 
in terms of n, which has one real solution 

V(r„) 2(3i 
n — 



3« V(r„)' 
where V(r) is given by 

V(r) = 3^/3 {^r^a'i + ^/?>^27r''a^e + ^a^P^eS^'^ 
We finally get the expression for potential 

^ ^" ~ V(0 ~ 2{V{rf-QaKf 
If we put /3 = in (j31|) the potential takes simple form 



Substituting (jH^ to (fT^ - (fT3j) we get correct energy levels (J2H1)- It is inter- 
esting to provide the expression (j^ when all constants are substituted, we 
have 

^^, , 3 / huTTy/m^ 
U[rj = -\r 



2 V r 

6 General Case of Arbitrary Spectrum 

One can find explicit relations for the radii of the orbits for arbitrary spectrum 
En. Indeed, From (J14j) and (J15p we want to find f/ as a function of r in terms 
of a given energy spectrum En. As we did above we assume the continuity 
of the parameter n and take derivative of both parts of (fT3j) in n, we have 

n n^2 n^ ^ "^ dr 



r=rn 



where prime denotes the derivative in n. The use of (fT^ allows us to get a 
differential equation for r„ in terms of only known quantities. We have 

/ 1 2 '^^ 771/ 

n 4 

11 



Introducing a new variable p = r"^ we obtain a linear differential equation 
which could be rewritten as 

-n L„ 



dn \ n / ^ 
which can be integrated to obtain 




n 



nEn -El- Ekdk + e (33) 



Equation ()33|1 expresses n-th radius in terms of n, i.e. it has the form r = 
f{n) now if we invert it we relate n in terms of r, n = /~^(r), which if 
substituted to (fT3j) gives an equation for U in terms of r only (actually in 
terms of r„, but we perform interpolation effectively ignoring the fact that 
the relation strictly holds only for orbit radii). 

Note that in (|33p we have a constant e (having units of energy) arising 
due to the integration, this means that we have a set of potentials resulting 
in the same energy spectrum. As we see from (jH^ the constant e could be 
determined if for example the smallest radius ri is known. In sections lU and 
El the situation was the same resulting in the constant P (see (E!?|l and (jSI})- 
Since e is more interesting from the point of view of physical interpretation 
we provide the expressions for (3 in terms of e. For the case of energy levels 
of hydrogen atom we have (see 



f^hydr = ^ ln(e + 7) 
and for the case of harmonic oscillator we have (see ()31j] ) 

Pose Z 



7 Energy Spectrum in Various Potentials 

As we already seen, for a given potential it is straightforward to compute 
corresponding energy levels. Indeed, from ()14j) we find r = Vn and upon 
substitution to ()15|l we get En- Like in quantum mechanics, potential U = 
U{r) should be attractive and strong enough to result in finite motion. Below 

12 



we consider several common central potentials which result in rather simple 
expressions for energy spectrum. In what follows the constant ^ is given by 
(fTH|) . note that it depends on time discreteness parameter as 1/r^. 
a). Coulomb potential 

Note that energy spectrum is different from the — 7/n^ spectrum obtained 
in quantum mechanics for this potential (see section on energy levels of hy- 
drogen atom for detailed discussion). 

b). Linear potential 

U[r) = ar, Vn = ——, E„ 



e ' " 2^ 



c). Logarithmic potential 



U{r) = alnr, rn = nJ-, E^ = a 



\ + '" («Vf 



For potentials (a), (b), and (c) r„ > if a > in all three cases it corresponds 
to the attraction field. 

d). Polynomial potential 

1 n- 

f/(r) = ar^ r„ = \-^\ , ^n = 2«(2 + ^) f ^ ) (34) 

This case generalizes cases (a) and (b) described above, although because of 
importance of these potentials we provided corresponding expressions explic- 
itly. Note that if we make a in (jMj) satisfy the following equation 

^^ 1, 



2-0- 



i.e. o = 2/3 then we get the linear energy spectrum En ^ n for the 2D 
quantum harmonic oscillator (see previous section for detailed discussion). 
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8 Discussion and Conclusion 

In this paper we have shown that discrete time formahsm leads to some 
distinguishable properties of micro-observables that are used to be described 
with quantum mechanics. In particular it was shown that finite motion 
results in discrete energy spectrum. Of the main interest in this paper is 
discrete energy levels of hydrogen atom. We have shown that for unperturbed 
hydrogen atom the discrete time formalism is able to give correct energy 
spectrum, more precisely we have reconstructed the corresponding "micro" - 
potential. Here we did not consider Stark or Zeeman effects for hydrogen 
atom, it would be interesting to study them from the point of view of discrete 
time formalism. 

As we have seen in above the discrete time model requires potentials which 
are different from QM potentials. One may argue this as a disadvantage of the 
model. We pay attention that there are no reasons to expect to reproduce QM 
by using the standard classical potentials; Bohr, Zommerfeld, Heisenberg and 
many others tried to do this, but they did not succeed. D. Bohm developed 
a new model [201 i^ which quantum mechanics can be reproduced on the 
classical basis, but, of course, the classical potential could not be preserved - 
it is perturbed by the quantum potential. And the latter looks not so natural 
from the classcial viewpoint, see e.g. the quantum potential for the two slit 
experiment in [201 ■ The discrete time model has an analogy with Bohmian 
mechanics - it tries to reproduce QM by changing potentials. But there is 
of course the fundamental difference: the only postulate that is used in the 
proposed approach is that there exists a quant of time r. 

Another interesting point is that one might expect that our dynamical 
equations are essentially difference equations which might produce discrete 
spectrum. This is not correct - we recall that in our model only time is 
discrete, but space is still continuous. 

There is still an open problem of the quantative value of the discreteness 
parameter r. One might speculate its relation with Plank time constant 
[S3 inSl ^ the smallest measurable time interval in ordinary QM and gravity 
- which is quantatively given by 



hG 

I' 



tpi = J'^^ 5.3910"^^ (sec.) 



The detailed analysis of this issue is out of the scope of the present article. 



14 



There might be an interesting interconnection on how the discrete time is 
used in information dynamics theory [TIH IT^ and the discrete time dynamics 
as it appears in presented study. In particular, it would be interesting to 
consider the equations for information dynamics with discrete time. 

Finally, we would like to comment also that there might be a deep in- 
terrelation between the energy-time uncertainty relations |28| 1221 ISHI and 
Bohr-Somerfeld quantization rules |22] in quantum mechanics and our dis- 
crete time model. In particular one can try to write the Bohr-Somerfeld 
semi-classical quantization rules for energy and time as canonical variables. 
For a system with conserved energy one might get EnTn ~ nH, this relation 
holds for example for energy levels En and classical periods T„ of Hydrogen 
atom. One the other hand the relation (|12p in discrete mechanics could be 
treated as the condition for the period to take only discrete values T„ ~ nr. 
We can see that although relations are similar there is an extra factor En in 
the quantum-mechanical relation. In fact one may argue that if we make the 
T in equations of motion depend on the energy of the system as 

e 

where tq is the "fundamental" time quantum and e a "fundamental" energy 
quantum, we get precisely the semiclassical quantization rules. The question 
arise how to treat the energy E here and what will happen with the dynamics. 
Further investigation of this interrelation will be discussed elsewhere. 
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